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Abstract 

We study the process of the association production of chargino and neu- 
tralino including the NLO QCD and the complete one-loop electroweak correc- 
tions in the framework of the minimal supersymmetric standard model(MSSM) 
at the Fermilab Tevatron and the CERN Large Hadron Collider (LHC). In 
both the NLO QCD and one-loop electroweak calculations we apply the algo- 
rithm of the phase-space slicing(PSS) method. We find that the NLO QCD 
corrections generally increase the Born cross sections, while the electroweak 
relative corrections decrease the Born cross section in most of the chosen 
parameter space. The NLO QCD and electroweak relative corrections typ- 
ically have the values of about 32% and —8% at the Tevatron, and about 
42% and —6% at the LHC respectively. The results show that both the 
NLO QCD and the complete one-loop electroweak corrections to the pro- 
cesses pp/pp X1X2 + generally significant and should be taken into 
consideration in precision experimental analysis. 
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I Introduction 



People believe that the minimal supersymmetric standard model(MSSM) is a very 
attractive extension of the standard model(SM). The direct discovery of the super- 
symmetric particles is one of the most important endeavors of present and future 
high energy experiments. The MSSM theory predicts many new particles, such as: 
squarks, sleptons, neutral CP-even Higgs-bosons h^, and CP-odd Higgs boson 
A", charged Higgs-bosons H^, charginos = 1,2), which are the mixtures of 
charged winos and charged higgsinos, and four neutralinos = 1 — 4) being the 
mixtures of the neutral wino, bino and two neutral higgsinos. 

In most of the constrained MSSM scenarios, such as the R-conserving mini- 
mal supergravity(mSUGRA), the lighter chargino(x]'^) and the neutralinos (x?, X2) 
are considerably less massive than the gluinos and squarks over most of the pa- 
rameter space, they may belong to the class of the relative light supersymmetric 
particles. The association production of xf c'-iid X2, which seems to be one of the 
primary source of the trilepton events, is a promising channel for supersymmet- 
ric particle searches at hadron colliders. In the association production of the lighter 
chargino(xf ) and the second lightest nuetralino(x2), both final supersymmetric par- 
ticles could dominantly decay to leptons, i.e., xt ~^ Xi'^'^^i and X2 ~^ Xi^~^^~ ^ which 
leads to a gold-plated £'^£~^i~ trilepton signature pi P|. We expect that the lighter 
charginos and neutralinos to be detected in their above mentioned decays at the 
Tevatron and the LHC. Actually, such trilepton signature was exploited primarily 
in the CDF and DO experiments at the Tevatron and got the similar bounds to those 
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from LEP2 in the MSSM parameter spacejH]. 

To investigate the discovering potential of hadron colhders, not only a proper 
understanding of the hadron production mechanisms is necessary, but also accuracy 
theoretical predictions of the signature should be provided. As we know, the impact 
of higher order electroweak and QCD contributions normally grows with increasing 
colliding energy and would become to be more obvious at very high energies. In 
Ref. V. Barger and Chung Kao investigated the prospects for detecting trilepton 
events {I = e or /i) from neutralino-chargino(x2xf ) associated production at the 
upgraded Tevatron in the mSUGRA model. They found that if there is a large 
integrated luminosity(for example C = 30 fb~^) and the decay modes Xi ~^ Xi^^^i 
and X2 ~^ Xi^^^~ kinematically dominant, the value of statistical significance 
(Ns = S/ ^/B, S =number of signal events, and B =number of background events) 
can reach 36.9 when we use suitable cuts and take tan/? = 3. Then one can expect 
the accuracy of the cross section measurement of the neutralino-chargino(x2Xi') 
associated production at the upgraded Tevatron can reach few percent. Therefore, 
for the precise experiments at TeV scale hadron colliders, both the higher order QCD 
and electroweak corrections should be considered in the theoretical predictions, and 
thereby one can improve experimental mass bounds and exclusion limits for the new 
particles. Moreover, the consideration of higher order QCD contributions can reduce 
the dependence of the cross sections on the renormalization and factorization scales 
in the LO. And the cross sections in NLO are under much better theoretical control 
than the leading order estimates. 

There have been many works which present the theoretical calculations of the 
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production of SUSY particles in hadron collisions at NLO QCD level, such as, 
Refs. |3] P [7| |H] . In Ref.jHj it presents the complete next-to-leading order SUSY 
QCD analysis for the production of all possible pairs of noncolored supersymmetric 
particles, including pp/pp X1X2 + their calculations of the NLO QCD 

corrections, the infrared and coUinear singularities were extracted by applying the 
dipole subtraction method jH]. 

In our work, we are to calculate and discuss the complete one-loop electroweak 
radiative corrections to the processes pp/pp xfx2 + at the Tevatron and the 
LHC. And for the completeness of our investigation we present also the NLO QCD 
corrections to the pp/pp xfx2 + processes by applying the phase-space slicing 
methodfPSS)[TU]. At the same time we compare our NLO QCD numerical results 
from the PSS method with those from dipole substraction method in Ref.jH]. 

The structure of this paper is organized as follows: In Section 2, we calculate 
the cross section of the leading order results for the subprocess ud — > XiX'2 ■ 
Section 3, we give the analytical and numerical calculations and discussions for the 
NLO QCD corrections to the processes pp/pp — > XaX2 + at the Tevatron and the 
LHC. In Section 4, we present the calculations of the complete one-loop electroweak 
corrections and discussions for the process pp/pp — > xfx2 + ^- Finally, a short 
summary is given. 
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II The leading order calculation for subprocess 



Since the cross sections for the subprocess ud — > X1X2 its charge-conjugate 
subprocess ud —>■ X1X2 CP-conserved MSSM are the same, we present here 

only the calculation of the subprocess ud — > X1X2 ■ The tree-level diagrams for the 
subprocess ud — > X1X2 ^^re shown in FigHJ 



Figure 1: The tree-level Feynman diagrams for the subprocess ud — > Xi X2 ■ 

In our calculation, we neglect the small masses of the light-quarks and there is 
no contributions from the Feynman diagrams which involve the couplings /G^ — 
u — d. 

We denote the subprocess as 



where Pi {i = 1, 2) and ki {i = 3, 4) are the four-momenta of incoming up-quark/anti- 
down-quark and outgoing xt and X2, respectively. All these four-momenta satisfy 
the on-shell conditions: /c| = m'i+, k\ = m~o, Pi = pi = 0- The center-of-mass 



ud 





u{pi) + d{p2) ^ Xtih) + xlih) 



(2.1) 



energy squared is denoted by s = E"^^ 



cm 



{Pl+P2f- 
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The amplitude of ud X1X2 subprocess can be divided into three parts and 
expressed as 

A^o = M' + M' + M"" (2.2) 

where Ai^ is the tree-level amplitude and Ai^, Ai^, M.^ represent the amplitude 
parts arising from the s-channel, t-channel and u-channel diagrams shown in Fig^a- 
c), respectively. Then the lowest order cross section for the subprocess ud —>■ X1X2 
in the MSSM is obtained by using the following formula: 

a,{s,ud^ Xtxl) = TTT-^ / diT\M'\'. (2.3) 

where 

imin{imax) = ^ { ("^~+ + "^|o ' s) ± ^ {s - m'i+ - m|o)2 - Am'i+m'ia | . (2.4) 

The summation is taken over the spins and colors of initial and final states, and 
the bar over the summation denotes averaging over the spins and colors of initial 
partons. 

Ill NLO QCD corrections 

In considering the NLO QCD corrections to the subprocesses ud{ud) XaX^-, we 
should involve the gluon emission subprocess ud X1X2 +fi' to cancel the soft IR 
divergence arising from the virtual QCD corrections of the subprocesses ud{ud) —>■ 
xfx2- The light-quark emission subprocesses gu{d) xtX2 + ^-i^) should be also 
included for a consistent and complete mass factorization. 
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III.l Virtual corrections 

The complete one-loop Feynman diagrams of QCD corrections in the MSSM, which 
are built up with gluon, gluino, quark and squark exchanging loops, are shown in 
Fig 121 We use the fermion flow prescription^^ for the calculation of the matrix 
elements including Majorana particles. Same as in the tree-level calculation, we 
neglect the masses of the light-quarks in the calculation of the virtual QCD correc- 
tions for the subprocess ud ^ X1X2 ■ There exist both ultraviolet (UV) divergences 
and soft /collinear IR singularities in the amplitudes from these QCD one- loop dia- 
grams. We use the method of the phase-space slicing(PSS) jTU] to treat the soft and 
collinear divergences. The PSS method is intuitive, simple to implement, and relies 
on a minimum of process dependent information and has been used in many works. 
In the NLO QCD calculation, we adopt the 't Hooft-Feynman gauge and dimen- 
sional regularization method with n = 4 — 2e to evaluate the one-loop contributions. 
The expressions for the relevant renormalization constants can be found in next 
section, except all the relevant self-energies including only the QCD parts instead 
of the electroweak parts. We have verified the cancellation of the UV divergence in 
the virtual QCD corrections analytically. Then we get an UV finite amplitude for 
the 0{as) virtual radiative corrections. 

The virtual correction to the cross section can be written as 



where is the UV renormalized amplitude for virtual corrections, and again the 

7 




(3.1) 
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summation with bar over head means the same operation as before. The expressions 
for imax and imin can be found in Eq. ()2.4|) 

After the renormahzation procedure, is UV-finite. However, it still contains 
the soft /coUinear IR singularities. The IR divergence part in da^ can be obtained 

as 



The soft divergences can be cancelled by adding the soft real gluon emission 
corrections. The coUinear divergences together with those coming from the real 
light-quark emission corrections are absorbed into the parton distribution functions, 
which will be discussed in the next subsection. 

III. 2 Real emission corrections 

The real emission subprocesses, which present NLO QCD corrections to the ud — >■ 
subprocess, include real gluon emission subprocess ud X1X2 +5'? and real 
light-quark emission subprocesses gu X1X2 + '^5 gd ^ X1X2 + ^- The later 
two subprocesses will make additional contributions to the final two-body and the 
three-body cross sections. 

1. Real gluon emission corrections 

The real gluon emission subprocess ud — > X1X2 presents 0{as) corrections 
to the subprocess ud — ^ X1X2 ■ is also one of the origins of IR singularities. Its 




(3.2) 



where 




= -4. 



(3.3) 



9 



IR singularities can be either of soft or collinear nature and can be conveniently 
isolated by slicing the ud — > X1X2 phase space into different regions defined 
by suitable cutoffs, a method which goes under the general name of phase-space 
slicing(PSS)flO . The soft IR singularity part from subprocess ud X1X2 +5' 
cancels exactly the analogous singularity presented in the 0{as) virtual corrections 
calculated in above subsection. 




(e) ^ (f) 

Figure 3: The tree- level Feynman diagrams for the real gluon emission subprocess 
ud xt')& +g- 

We denote the 2^3 subprocess ud — > X1X2 +5' 

u{pi) + d{p2) ^ xtih) + xlikA) + g{h), (3.4) 

The Mandelstam variables are defined as 

s = (pi + p2)^ i = {pi - ksY, u= {pi- kif, S45 = {ki + k^f, 
ii5 = (Pi - hf, ^5 = (P2 - k^y, £45 = {ki - k^f (3.5) 

By introducing an arbitrary small soft cutoff 5s we separate the phase space of the 
2 — > 3 subprocess ud X1X2 i'^^o ^"^o regions, according to whether the energy 
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of the emitted gluon is soft, i.e. < or hard, i.e. > The 

partonic cross section can be written as 

afiud^ Xtxl + 9) = oliud^ xtxl + 9) + ^f{ud-^xtxl + 9), (3-6) 

where is obtained by integrating over the soft gluon emission phase space region, 
and contains all the soft IR singularities. In order to isolate the remaining collinear 
singularities from , we further decompose into hard-collinear(HC) and hard 



non-collinear (HC) parts by introducing another cutoff 6c named collinear cutoff, 
^fi^d ^ xtxl + 9) = ^f'iud ^ xtxl + 9) + ^'^{ud xtxl + 9)- (3-7) 

The HC regions of the phase space are those in collinear condition, where any 
invariant, ti5,t25,^45, becomes smaller in magnitude than 5c V^, while at the same 
time the emitted gluon remains hard, a^'^ contains the collinear divergences. In 
the soft and HC regions, and a^^ can be obtained by performing the phase 



space integration in n = 4 — 2e dimensions analytically. But in the HC region, a^^ 
is finite and can be evaluated in four dimensions by using standard Monte Carlo 
techniques jl2j. The cross sections, d^, a^'~^ and o"^*", depend on the arbitrary 
parameters, 5s and 5c. 

With the arbitrary small cutoff 5^, the differential cross section in the soft region 
is given as 



da^g = da° 



with 



Af = | Af = -^ln5„ Ao^ = yln^5.. (3.9) 
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In the limit where two of the partons are colhnear, the three body phase space 
is greatly simplified. And at the same limit, the leading pole approximation of the 
matrix element is valid. According to whether the colhnear singularities are initial 
or final state in origin, where a^^ can be separated into two pieces 

*r=*,T+*,"f (3.10) 

Since in the process ud — > XiX% +§> only initial particles are involved in strong 
interaction and massless in the limit, we only need to calculate the cross section 
which arises from the case that the emitted gluon is colhnear to the initial 
partons, i.e., < ^15,^25 < ^cS. 

In HC phase space region, the initial state partons i{i = u,d) is considered 
to split into a hard parton i' and a coHinear gluon, i i'g, with pi' = zpi and 
fcs = (1 — z)pi. The cross section a^/" can be written as 

<^ = daff = da' ^^^^^^ i^^Jj (^-^^ 

dz fl-z\~'^ 

where G^ j/a,b the bare parton distribution functions. A and B refer to protons 
at the LHC, and proton, antiproton at the Tevatron. Pu„(z,e) and Pjj{z,e) are the 
n-dimensional unregulated {z < 1) splitting functions related to the usual Altarelli- 
Parisi splitting kernels [T^. Pii{z,e){i = u,d)ca.n be written explicitly as 



Pu{z,e) = Pu{z) + ePl,{z) 
Piiiz) = Cp 



1 + z' 



1 - z 

Pi{z) = -Cf{1-z) {z = u,d), (3.12) 
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with Cf=4/3. 

2. Real light-quark emission corrections of gu{d) X1X2 + d{u) 
In addition to the real gluon emission subprocess, the subprocesses gu{d) 
XiXH + ^{^) should also be included. The contributions from these two subprocesses 
contain only the initial state collinear singularity. These subprocesses will make 
additional contributions to the two-body term cross section o"'-^^ and the three-body 
term cross section a'-'^^ The Feynman diagrams for these two subprocesses at the 
tree-level are shown in Fig0]and FiglHl respectively. 




(e) (f) (g) (h) 

Figure 4: The tree-level Feynman diagrams for the real light-quark emission sub- 
process gu xtxl + d. 



By using the PSS method described above, we split the phase space into two 
regions: collinear region and non-collinear region. 

af{gu{d)-.xtxl + d{u)) = af''{gu{d)^xtx'2 + d{u)) 

+ af^{gu{d)^xtX2 + d{u)). (3.13) 
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(a) 



(b) 



(d) 




(e) (f) (g) (h) 

Figure 5: The tree-level Feynman diagrams for the real light-quark emission sub- 
process gd X]l"x2 + 



Also in hard non-coUinear region is finite and can be evaluated in four di- 



mensions using Monte Carlo method. The differential cross sections of da^ for 
subprocesses gu{d) Xi)(^ + d{u) can be written as 



daf'^{guid)^xt)^, + diu))^ 



da' 



27rr(] 

+ {A^B)]—[- — -] dxAdxB 



\ ) ^c'[Pd{u)g{z,e)Gg/A{xA/z)G^^d)/B{xB) 



z \ z 



(3.14) 



with Pij{z, e){i = u, d, j = g) expressed explicitly as 



Pd{u)9{z,e) = Pd(u)g{z) + eP'^^^-^g{z), 
PdiuUz) = liz' + il-zf], P^(^)^(^) = -^(l-^). (3.15) 

3. NLO QCD corrected cross section for pp/pp + ^ 

After adding the renormalized virtual corrections and the real corrections, the 
partonic cross sections still contain the coUinear divergences, which can be absorbed 
into the redefinition of the distribution functions at NLO. Using the MS scheme, 
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the scale dependent NLO parton distribution functions are given as 



as r(l — e) / 47r/i, 



27rr(l-2e) I /i^ 



^ dz 



Pij{z)G.j/A{x/z). 



(3.16) 



By using above definition, we get a NLO QCD parton distribution function counter- 
terms which are combined with the hard colhnear contributionsfEq . (j3. 1 Ij) and (j3.14j) ) 
to result in the 0{as) expression for the remaining collinear contributions: 

coll J^O \(^S -e) ( 4:71 fij^ 



da 



da" 
E 



27rr(l -2e) V s 
Af{a — > ag) 



{Gu/a{,Xa, IJ,f)G^/BixB, fif) + Gu/AixA, IJ,f)Ga/BixB, yU/) 



+ A'o'ia ag) 



a=u,d 

{A ^ B)}dxAdxB, 



Gu/a{xa, fIf)G^/B{xB, fJ'f) 



where A and B are proton and antiproton for the Tevatron, and proton and proton 
for the LHC, respectively. 

Al'iuid) ^ u{d)g) = GF{2\n6s + 3/2), A'^' = Al" In I -^V (3.18) 



and 



Ga/A,B{^^f^f) = 



c'=a,g 



dy 



Gc'/A,B{x/y,fif)Pac'{y), {a = u,d) (3.19) 



with 



PaAy) = Pac In I 5c^4 1 - KAy)- 

y f^f 



(3.20) 



We can observe that the sum of the soft fEq. ()3.8p ). collinear ( Eg . ()3 . 1 7|) ) . and 
ultraviolet renormalized virtual correction (Eq. ()3.2|) ) terms is finite, i.e.. 



Al + A^ = 0, 



Af + A^ + Af{u ug) + Af{d dg) = 0. 



(3.21) 



(3.17) 
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The final result for the total 0{as) correction consists of two parts of contributions: 
a two-body term cr*-^-' and a three-body term a^^\ The two-body correction term 
0"^^) is expressed as 

^^^^ = ^ J dxAdxBda^{Gu/AixA,l^f)Ga/B{xB,f^f)[A^ + A'^ + A'Q''{u-^ug) + A'o^{d^dg)] 
+ Gu/AixA, /i/) + Gu/AixA, lJ-f)G^/j^{xB, fif) + (A ^ B)}. (3.22) 

And the three-body correction term is written as 

+ (T^^) (pp/pp ^ gu^ xtxl + d) + (t(^) (pp/pp ^ gd^ xt^2 + "«) 
= y dxAC?a;B[G„/^(xA,/i/)G'j/B(a;B,/i/) + (A ^ 5)]rf(T(^)(nJ^ Xi^x^ + 5') 
+ j dxAdxB[Gg/A{xA, f^f)Gu/BixB, fJ'f) + {A ^ B)]da^^\gu xtxl + d) 
+ J dxAdxB[Gg/A{xA,fJ^f)G^/B{xB,fJ'f) + {A^ B)]da^^\gd^xtx^2 + u). 

(3.23) 

Finally, the NLO total cross section for pp/pp X1X2 + is 

In the cross section part of a^'^^ + cr^^\ the dependence on the arbitrary cutoffs Sg 
and 6c should be vanished. This constitutes an important check of our calculation. 

III. 3 Numerical results involving NLO QCD corrections 

In Ref.|Hj, W. Beenakker, et al, presented the NLO QCD calculations of the processes 
pp/pp — » X1X2 + There the infrared and coUinear singularities of the three 
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parton cross sections are extracted by applying algorithm of the dipole subtraction 
methodjn]. As a check of our numerical calculation by adopting the two cutoff 
PSS method, we take the same mSUGRA input parameters [mi/2 = 150 GeV, uiq = 
100 GeV, Ao = 300 GeV, /i > 0, tan /? = 4] and reproduce the Fig.2 shown in Ref. 
with the coincident numerical results. From this comparison we have verified the 
correctness of our calculations of the NLO QCD corrections. We take the colliding 
energies of pp/pp at the Tevatron Run 11 and the LHC are 2 TeV and 14 TeV, 
respectively. 

In the following numerical calculation for NLO QCD corrections, we use the 
package FormCalc J3] to get all the masses of the supersymmetric particles by in- 
putting tan/3, rriAo, Mgusy, A*, ^2 and Af parameters. In the package FormCalc, 
the grand unification theory(GUT) relation Mi = (5/3) tan^ 6'vi/M2 is adopted [T^. 
and Mq = Mjj = Mjj = = Mi = Mg^sy is assumed in the sfermion sector for 
simplification. We use the one-loop formula for the running strong coupling constant 

with as{mz) = 0.1187, and take CTEQ6L and the CTEQ6M parton distribution 
functions in calculating the LO and the NLO cross sections, respectively^?]. 

Figini shows the independence of the NLO QCD corrected cross sections on the 
arbitrary cutoffs Sg and 6c by applying the two cutoff PSS method. The two-body 
correction term (cr*^^^) and three-body correction term(o"*^'^-*) and the NLO QCD 
corrected total cross section a'^^^ at the Tevatron and the LHC, are shown as 
the functions of the soft cutoff 6s with the collinear cutoff 6c = 6s/ 50 (shown in 
Fig|nfa),(c)), and as the functions of the collinear cutoff 6c with the soft cutoff 
6s = 506c (shown in Fig|nfb),(d)). The input supersymmetric parameters are taken 
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as tan/3 = 4, m^o = 300 GeV, M^usy = 250 GeV, /i = 278 GeV, Ma = 123 
GeV and Aj = 450 GeV. We can see the NLO QCD corrected total cross section 
(jNLO _ ^o_j_^^ _ ^o_|_^(2) _|_^(3) jg independent of the cutoffs. This is an important 

check of the correctness of our calculation. In the following numerical calculations, 
we set 5s = 10~^ and 6c = 6s/ 50. 

During our numerical calculation we investigate also the dependence of the LO 
and NLO QCD corrected total cross sections at the LHC and the Tevatron on the 
renormalization and factorization scales(Q = fir = f^f), and find that the theoretical 
predictions including NLO QCD corrections become stable, being nearly indepen- 
dent of the factorization/renormalization scales for the processes pp/pp xfx2~^^ 
as concluded in Ref.jHJ. 

In FigQ the dependence of the Born, NLO QCD corrected cross sections (shown 
in Fig|7|^a) and (d)) and the relative corrections 6 (shown in Fig|7U^b) and (d)) for 
the processes pp/pp — > xfx2 + at the LHC and the Tevatron on the gaugino 
mass parameter M2 are plotted. There we take the input parameters as m^o=300 
GeV, M,„,j,=500 GeV, /i=400 GeV and Af=m GeV, with tan /5 = 4, tan /3 = 15 
and tan/5 = 40, respectively. We can see the Born and NLO QCD corrected cross 
sections decrease rapidly to a small value with the increment of M2. At the Tevatron, 
the NLO QCD relative correction for tan/5 = 4 decreases from 37.7% to 29.0% with 
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Figure 6: The dependence of the total cross sections forpp/pp — > XiX2^^ processes 
at the Tevatron and the LHC as the functions of the cutoff 6s with 6c = 6s/ 50 (see 
Fig\^a),(c)) and 6c with 6s = 506c(see Fig\^h),(d)), respectively. 



19 



I 1 1 1 1 1 
(a) 


§LHC , 
















LO 




NLO 


I.I.I 





120 



120 



140 160 

M(GeV) 



1 1 1 1 
@LHC (b) 


1 1 ' 


^. 








'\ tg(J=4 

V \ / 




■\ X 
\ \ 




\ \ 




\ \ 
\. \ 




\ \ 
'v \ 
s,. \ 


tgp=40 


s- \. 
s ■ 


z : 


tgp=15 

I.I.I. 



180 



200 120 



140 



160 

M,(GeV) 



180 



200 



1 1 1 1 1 
(c) 


1 1 
@Tevatron- 


■■. ,tgp=4 






V.', '-y 






\^tgp=15 
















LO 






NLO' 


1 1 1 1 1 


1 





^ 34 - 



©Tevatron 


1 1 1 1 1 
(d) 




\tap=4 








tgp=46-.^^^C\ 


■ 


I.I. 



140 160 



120 



140 



160 



180 



200 



180 200 

(GeV) (GeV) 

Figure 7: The dependence of the Born, NLO QCD corrected cross sections (shown 
in Fig. 12(a)) and the corresponding relative corrections 5 (shown in Fig. 12(h)) for 
the processes pp / pp xfx2 + ^ '^^ LHC and the Tevatron on the gaugino mass 
parameter M2. There we take the input parameters as mA0=300 GeV, Msusy=500 
GeV, 11=400 GeV and Af=450 GeV, with tan/3 = 4, tan/3 = 15 and tan^S = 40 
respectively. 
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the increment of Ma from 120 GeV to 200 GeV. While at the LHC, the relative 
corrections decrease at the region M2 < 160 GeV, and then go down slowly in the 
region M2 > 160 GeV. The NLO QCD relative corrections for tan/5 = 4, 15,40 at 
the LHC have the values in the range between 45.3% to 41%, when M2 runs from 
120 GeV to 200 GeV. 



In FiglHfa) and (c), we present the Born cross sections and NLO QCD corrected 
cross sections for the processes pp/pp xf X2 + X aX the Tevatron and the LHC 
as the functions of the supersymmetric soft breaking mass parameter Mgusy, on the 
conditions of m^o=300 GeV, M2=123 GeV, ;U=278 GeV and A/=450 GeV, with 
tan/5 = 5, tan/5 = 20 and tan/5 = 40, respectively. The solid and dotted curves 
present the Born and NLO QCD corrected cross sections, respectively. We can 
see the cross sections, especially the corrected ones, increase with the increment of 
Msusy FiglHJb) and (d) present the relative corrections 5 as the functions of Msusy 
at the LHC and the Tevatron, respectively. The relative correction at the Tevatron 
for tan/? = 4 increases rapidly from about 29.7% to 38.4% when Mg^sy goes from 
250 GeV to 600 GeV , and becomes to be a constant value about 38.4% when Mgusy 
varies in the region beyond 600 GeV. While at the LHC, the relative correction 
decreases rapidly when Mg^sy goes up from 250 GeV to 400 GeV , and keeps the 
value about 45% in the region of Msusy > 400 GeV. 

We also calculate the cross sections of the processes pp/pp xfx2 + as the 
functions of higgsino-mass parameter /i(with m^o = 300 GeV, M2 = 120 GeV, 
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Figure 8: The Born cross sections, NLO QCD corrected cross sections and the 
relative NLO QCD corrections S for the processes pp/pp xfx'i'^^ o,t the Tevatron 
and the LHC as the functions of the supersymmetric mass parameter Mgusy, on the 
conditions ofmAo=300 GeV, M2=123 GeV, ^=278 GeV and Af^450 GeV, with 
tan P — 4, tan P — 20 and tan P — 40, respectively. 



22 



Msusy = 500 GeV, Af = 450 GeV, tan/3 = 4 or 15 and /i E [250 GeV, 1000 Gel/]), 
and tan/?(with m^o = 300 Gel/, M2 = 160 GeV, Msusy = 500 GeV", Af = 450 GeV, 
fi = 450 Gel/ and tan/5 G [5, 40]) at the LHC/Tevatron and find the result does 
not depend much on these parameters. And the relative corrections have the typical 
values of about 32% and 42% at the Tevatron and the LHC, respectively. 



IV Full one-loop electroweak corrections to the 
ud X1X2 subprocess 



For a precise analysis of the association production processes of neutralino and 
chargino, higher order electroweak corrections should be included. The one-loop 
level UV renormalized electroweak virtual corrections to the ud -—>■ X1X2 subprocess 
can be expressed in the form as 



And again the summation with bar in the equation recalls the same operations as 
appeared in Eq. ()2.3|l . imax and imin are the same as expressed in Eq. ()2.4j) . Ai^ 
is the the UV renormalized amplitude of virtual Feynman diagrams including self- 
energy, vertex, box and counterterm diagrams. We use FeynArts and FormCalcfT^ 
packages to generate Feynman diagrams. For the numerical evaluation of the loop 
integrals we use the developed package LoopTools^^. 

IV. 1 Renormalization scheme 

As we know, the contributions of the electroweak one-loop diagrams contain both 
ultraviolet (UV) and infrared(IR) divergences. The UV divergence can be regu- 




(4.1) 
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larized by adopting the dimensional reduction(Di?) regularization scheme PU] and 
the relevant fields are renormalized by using the on-shell(OS) conditions |21j |26j 
(neglecting the finite widths of particles). In treating the QED soft and collinear 
divergences we use again the two cutoff PSS methodlTU], the IR divergencies are 
cancelled in complete analogy to the calculation for the corresponding QCD radia- 
tive corrections as shown in Section 3. The amplitudes are performed by adopting 
the 't Hooft-Feynman gauge and n = 4 — 2e space-time dimensions to isolate the UV 
and IR singularities. There is no QED induced collinear IR singularity from final 
state radiation due to chargino being massive, but there exists collinear IR singular- 
ities from the initial state radiation. Similar with that as declared in Section 3, in 
the electroweak correction calculation we use again the fermion fiow prescription to 
deal with the matrix elements including Majorana particles |TI]. The quark mixing 
matrix is assumed to be diagonal. The bare parameters are split into renormalized 
parameters and their counter terms. 
1. Gauge sector 

The definitions and the explicit expressions of the renormalization constants for 
the gauge boson sector are written as |^ 



^{l + ]^5Zw)W^, (4.3) 

Z\ (l + \6Zzz \6Zzz\(Z\ 
A ) \ \5Zaz 1 + )\a)- ^^■^> 

The renormalization conditions are taken as that the renormalized mass param- 
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eters of the physical particles are fixed by the requirement that they are equal to the 
physical masses, i.e., to the real parts of the poles of the corresponding propagators. 
Then the relevant fields and mass parameters are properly normalized. It yields the 
following results for their counter terms. 



6ml = Re^T^iml), bm\^ = ReE^^ {mly) , (4.5) 

Qp2 



SZvv = -Re 2 |p^=o, V = A,Z, W, (4.6) 



2i?eS^^(m|) 2i?eS^^(0) 

dZAZ = o 5 oZzA = 2 l^-') 

where denotes the transverse self-energy and Re means only taking the real part 
of the loop integrals. The weak mixing angle is define as = 1 — m^/m||2H]- 
2. Fermion sector 

The relevant fermion field renormalization constants are defined as: 

mf,o = mf + 6mf, f^^ = {I + Uz^^,^)f\ = {I + Uzl^) . (4.8) 

For the SM fermions, the normalized constants can be expressed as: 

5mf = ^Re m/Sj'(mJ) + m/S^(mJ) + Sj'^(mJ) + ^^'^(mj) , (4.9) 



6Z^,, = -ReT.l,,{m)) - mj—Re {m^Sj,,,(p2) + m^S^,,(/) 



d_ 

sJ,i(/) + sJ,;?(p^)} (4.11) 
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TTif — m. 



-Re 



R 



Ji Jj 



^ j] 



(4.12) 



m"^ — m. 



-Re 



Ji Jj 



K. ) + ; 



(4.13) 



The one-particle irreducible two-point function iV f^ij{p'^) for fermions is decomposed 



as 



T/., (/) = 5.. (l^ - rrif) + [i^P^E^,,. (p^) + ^P^E« (p^) 



+ ^Sj^(p') + ^Sj;5(p^ 



(4.14) 



In our calculation we use an improved scheme to make the perturbative calcu- 



lation more reliable. That means we use the effective MS fine structure constant 
value at Q = mz as input parameter, o.ew{jnz)~^\-MS ~ 127.918[34J. This results in 
the counter-term of the electric charge expressed as 1231201 EI]: 



Aj2Nle}iA + \og^]+J2J2Nle}iA + log 



f 



f 



f 



k=l 



2 

Jk , 



-22 A + lo£ 



m 



w 



(4.15) 



where we take xj = niz when mj < mz and Xt = rrit. e/ is the electric charge of 
(s)fermion and A = 2/e — •y + log47r. A*"^ is color factor, which equal to 1 and 3 for 
(s)leptons and (s)quarks, respectively. In our calculation we take Q = (m-±-|-m^o)/2 
in using Eq. ()4.15|) . 
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3. Supersymmetric sector 

In the MSSM theory the physical chargino mass eigenstates Xi^2 combi- 
nations of charged gauginos and higgsinos. Their physical masses can be obtained 
by diagonalizing the corresponding mass matrix X, which has the form as j37j : 

X=( ^su{2) \/2mw sin (3 \ 
\ \/2mw cos (3 J ' 

X is diagonalized with two unitary matrices U and V according to 

U*XV^ = diagim^^m^^), (4.16) 

which yields the chargino masses. 

The neutralinos are the mixtures of the neutral gauginos and higgsinos. Their 
physical masses can be obtained by diagonalizing the corresponding mass matrix Y 



Y 



/ Mu{i) —mzsw cos P mzSwsmP \ 

Msu{2) rnzcw cos (3 —mzcw sin (3 

—mzSw cos (3 mzCw cos P — /i 

\ mzSwsmP —mzCw sin P — J 



(4.17) 



Y is diagonalized with a unitary matrix N according to 

_/V*y_/Vt _ diag{m^o,m^o,m^o,m^o), (4-18) 

which yields the four neutralino mass eigenstates. 

We follow the renormalization definitions for chargino, neutralino and sfermion 
sectors as in Ref. |23j. The chargino, neutralino wave functions and mass counter 
terms in the mass eigenstate basis are introduced as 
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where x stands for both charginos and neutrahnos, i,j = 1,2 for chargino sector, 
i,j = 1,2,3,4 for neutrahno sector. These counter terms can be expressed as the 
functions of the corresponding self-energies similar with the equations of Eqs. (j4.9|) - 
()4.14|1 with the replacements of fij — ^ Xij. 

The wave function and mass counter terms for scalar quarks and scalar leptons 
are defined as 



1 + \m (I ) , (4.20) 



where 



SZ^^^^ I ^% 1. (4.22) 

\5Ztl 6Zg J 

The corresponding counter terms for the scalar fermions are given as 

Ha. = ^^"''(-L^.)' ^^t'' = (4.23) 

_2S|V_-) 

5Zll = -Re— ^- (^,j = l,2, 1^3). (4.24) 

ml - — mi - 

We introduce the counter terms for unitary matrices A^, U, V and R as follow: 

N^N + 6N, U^U + 6U, V^V + 6V, R^R + 6R, (4.25) 

where N and R are the rotation matrices of neutralino and squark(slepton) sectors, 
respectively. U and V are the diagonal unitary matrices for chargino sector. 
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The counterterms 5U, 5V, SN and SR can be fixed by requiring that the coun- 
terterms 6U, 6V, 6N and 6R cancel the antisjTumetric parts of the wave function 
correct ions Ell We get the expressions of the counter terms for the neu- 
tralino, chargino and sfermion rotation matrices A^, U, V and R as below: 

1 ^ ^ 

^^^^ = i - ^^kh^'^v (4-26) 

fc=i 

4 4 

^U^, = \ Y.{5Zt/ - bZt^)N,,, 5V,, = \ Y.{bZtf^ - 5^+'^)iV,, (4.27) 

k=\ k=l 

K = \ il^^zl - 6zI)rI (4.28) 

k=l 

As we expected, the UV divergence induced by the one-loop diagrams can be 
cancelled by that contributed by the counterterm diagrams exactly. While the soft 
and collinear IR divergences still exist. 

IV. 2 Real photon emission 

The soft IR singularity in the is originated from virtual photonic loop correction. 
It can be cancelled by the contribution of the real soft photon emission process. The 
real photon emission Feynman diagrams for the subprocess ud — > xtX2l shown 
in FigH 

We denote the real photon emission process as 

u{pi) + d{p2) xtih) + xUh) + 7(^:5), (4.29) 

We adopt the general PSS method|lO^ to separate the soft photon emission singu- 
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(e) (f) (g) (h) 

Figure 9: The real photon emission Feynman diagrams for the subprocess ud — >■ 

xtxli ■ 

larity from the real photon emission process. By introducing an arbitrary small soft 
cutoff 5s we separate the phase space of the subprocess ud X1X2I i^^o regions, 
according to whether the energy of the emitted photon is soft, i.e. < 6sV§/2, 
or hard, i.e. E^, > 6sV^/2, where is the incoming parton beam energy in the 

c.m.s. frame. Then the correction of the real photon emission is broken down into 
corresponding soft and hard terms 

A&real = soft + ^^hard = ^oi^soft + hard)- (4.30) 

Although both Ad soft and Aahard depend on the soft photon cutoff SsV^/2, the 
real correction Aareai is cutoff independent. If we use the soft photon emission 
approximation [221 ! we can set fc^ = in the delta function of the phase space 
element(i.e., 5^^\pi + P2 — — k^)) up to corrections of 0{Ss)- Then we can take 
the n-momenta of the initial and final particles in the pi + p2 rest frame as 

Pi = ^(i,...,o,o,i), p2 = ^(l,...,0,0,-l), 

^3 = (£"3, ...,psin6', 0,pcos6'), /C4 = (i?4, — psin^, 0, — pcos6'), 

/c5 = i?5(l, sin6'i sin 6^2, sin6'i cos6'2, cos6'i). (4-31) 

Then the integral over the soft photon phase space can be implemented analytically. 
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We get the expression of the differential cross section for the subprocess ud — > X1X2 
asdU 



dAa 



soft 



dcFQ 



OLe, 



r(i 



27r r(l -2e) 



(4.32) 



where 



A^ 
^0 



9' 

^ \ 12 ]^ 'o Y 

^2^*12 + "^^13 + gC'23 ^ C's: 



33' 



33' 



with 



12 



-2 In (5, 



^0 



^1 

'-'13 



^0 
'-'13 



^1 



'-^23 



'-'33 



'^33 



2W 5„ 

2 



In^ 



-E3 



^3 +P 

— p cos 6^ 2 — p 



+ 2Li 



p cos 6* — p 



^3 



- 2Li 



—p cos 6* — p 
E3 — p cos 6* 



(E3-pCOS^ 2 

+ Ind^ln -2 5 ^ ^^' 

1 (Es+pcos^)^ 



In^ 



^3-P 1 , 2 ^3 +P , or ~P ^^.PCOSO-P 

— In — \-2Li — 2L«— 

E3 + p cos 2 E3 — p E3 — p E3 + p cos 



{Es+pcosOy 2 . 
+ Ind^ln -2 r hln d^, 



21n5,-^ln.^^+^ 



-E3 



(4.33) 
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In above equations we have following relations between some variables: 



E3 



2m^+ - i - u 



p cos U = — 

2VI 




2VI 



p = 



2v^V 



sm-+ —i — u)i2m? 



2VIm-+ - i-u). (4.34) 



Since the incoming light-quarks are assumed to be massless in the parton model 
and the outgoing particles are massive, there exist only collinear IR singularities 
induced by initial state hard photon radiation. To isolate the collinear singularities 
from Adhard, we further decompose A&hard into a sum of hard collinear(HC) and 
hard non-collinear(HC) terms by introducing another cutoff 6c named collinear cutoff 



where the HC regions of the phase space are those any one of the Lorentz invari- 
ants ti5(= {pi — k^Y), %(= {P2 — k^Y) becomes smaller in magnitude than 5cS 
and the emitted photon remains hard. Acthc contains collinear divergences. As 
mentioned above, the soft IR divergence of the virtual photonic corrections can be 
cancelled exactly by that of soft real corrections. The remaining collinear singu- 
larities are absorbed by a redefinition (renormalization) of the parton distribution 
functions PDFs|j32j. This is done in analogy to the calculation of QCD radiative 
correction. 

In the HC region, Aa^^ is finite and can be evaluated in four- dimensions by 
applying standard Monte Carlo method. We can see that Adsoft, ^ctrc and Aa^, 
depend on the two arbitrary parameters Sg and 6c. However, in the total elecroweak 
corrected hadronic cross section, after mass factorization, the dependence on these 



(4.35) 
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arbitrary cutoffs(5s and Sc) cancels, as will be explicitly shown in numerical calcula- 
tion(see Section IV.4). This constitutes an important check of our calculation. Fi- 
nally, we get an UV and IR finite corrections Act of the processes pp/pp — > xf X2+^- 

IV. 3 The cross sections of processes pp/pp XiX^ + ^ 

The total electroweak corrected cross sections of the parent processes pp/pp 
ud{ud) — > xf + X at the Tevatron and the LHC can be calculated from the cross 
sections of subprocesses ud — >• x^X2 ud — >• Xi^- 

By summing the UV renormalized electroweak virtual corrections and the real 
photon emission corrections, the remaining coUinear divergences are absorbed into 
the redefinition of the distribution functions. Using the MS scheme, the scale de- 
pendent parton distribution functions including 0{aew) corrections are given as 



Gi/A{x,IIf) = Gi/A{x) + 



2n r(l - 2e) \ ~/if 



1 dz 

-P^,{z)G,iAx/z). 



(4.36) 



By using above definition, we get a 0{a.f.y^ parton distribution function counter- 
terms which are combined with the hard coUinear contributions to result in the 
0{af.y^ expression for the remaining coUinear contributions: 



+ E 



271 r(l -2ej \ s 
Af'{a — > 0:7) 



{Gu/a{xA: Hf)Gj/B{xB: l^f) + Gu/a{xa, Hf)Gj/B{xB: l^f) 



a=u,d 

+ {A ^ B)}dxAdxB, 



Gu/a{xa, i^f)Ga/B{xB, l^f) 



(4.37) 



where A/B are proton/antiproton for the Tevatron, and proton/proton for the LHC, 
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respectively. 

ATiuid) ^ uidh) = gJ(21n5, + 3/2), A^' = Ar\n[^], (4.38) 

and 

Ga/A,Bix,Hf) = 2^ / — qjGc'/A,Bix/y,Hf)Pac'{y), {a = u, d) (4.39) 
with 

Pac'(y) = ^ac'ln ( 5c^4 ) - ^ac'(2/)- (4.40) 

where Pqg = jPqq and Pg^ = ^Pqg, Pqq and Pqg have the expressions as shown in 
Eqs. fl3.12|l and ()3.15j) . respectively. g/(/ = u,d) are the charges of quarks. 

The final result for the total 0{aew) correction consists of two parts of contribu- 
tions: a two-body term cr*^^) and a three-body term a'^'^\ The two-body correction 
term cr^^'' is expressed as 

^(2) _ jdxAdxBda'^{Gu/A{xA,l^f)Gd/BixB,l^f)[A^ + A'^ + A'^''{u^u-f) + A'^''{d^d-f)] 

+ Gu/AixA, /i/)Gj/B(xB, /i/) + Gu/AixA, II f)G ^/ b{,Xb, /i/) + (A ^ B)} . (4.41) 

And the three-body correction term is written as 

^(3) ^ ^(3) (pp/pp ->ud--^ xtxl + 7) 

= j dxAdxB[Gu/AixA, ^f)Ga/B{xB, fJ'f) + (A ^ B)]da^^\ud xtxl + 1), 

(4.42) 

where G is the proton/antiproton distribution function. Finally, the full one-loop 
electroweak corrected cross section for pp/pp X1X2 + is 

= a° + + = a° + Aa. (4.43) 
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The cross section part of a^"^^ + a^^^ should be independence of the cutoff parame- 
ters 6s and 6c- The electroweak one-loop relative correction is defined as 6 = Aa/d-Q. 

TV A Numerical results including electroweak corrections 

In this subsection, we present some numerical results for the one-loop 0{aew) elec- 
troweak corrections to the processes pp/pp xf X2 + We take the SM in- 
put parameters as niz = 91.1876 GeV, my/ = 80.425 GeV, nit = 178.1 GeV, 
rrif, = 4.7 Gel^jSn] and neglect the light-quark masses in the numerical calculation. 
The fine structure constant is taken having the value at the Z°-pole, o:ew{fn'z)\j4s — 
1/127.918(311. The new MRST 2004-QED parton distribution functions including 
0{aew) corrections to the parton evolution are adopted in calculating the Born and 
the one-loop order corrected cross sections [35 . The renormalization and factoriza- 
tion scales are taken to be equal for simplicity((5 = /ir = A*/), and have the value 
being the average of the final particle masses in analogy to the NLO QCD calcula- 
tion. We use again the package FormCalc to obtain all the masses of supersymmetric 
particles by inputting the supersymmetric parameters tan/3, m^o, Msusy, fJ^, M2 and 
Af. Among these six input supersymmetric parameters, the CP-odd Higgs-boson 
mass m^o and tan/? with the constraint tan/3 > 2.5 are for the Higgs sector. In 
FormCalc package the radiative corrections to Higgs-boson masses up to two-loop 
contributions have been involved|3tij. While the tree-level Higgs-boson masses can 
be obtained by using the equations 




- ym\o + m|o =F ^J {m\o + "^|o)^ - 4m^o"^|o 
2 I 2 



cos2(2/5) 



(4.44) 
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In order to keep the gauge invariance during our numerical calculation, we adopt the 
tree- level Higgs-masses obtained from Eq. ()4.44|) , but not the Higgs- masses from the 
output of FormCalc package through out the tree-level and one-loop calculations. 

As mentioned above, the final results should be independent on cutoffs 6s and 6c- 
For demonstration, we present the cross section corrections of the pp/pp — > xfx2^^ 
processes as the functions of the soft cutoffs 6s and 6c in Figs llOr a-d) on conditions 
of tan/5 = 4, ruAo = 300 GeV, Msusy = 250 GeV, fi = 278 GeV, Ma = 127 GeV 
and Af = 450 GeV at the Tevatron and the LHC. The dashed, solid and dotted 
lines correspond to the total correction Aa = a^'^^ + a^^\ three-body correction 
a^^^ and two-body correction a^'^\ respectively. As shown in these figures, the full 
0{aw) correction Act is independent of the soft cutoff 6s{6c), as 6s{6c) running from 
10-6(10-'^) to 10-1(10-3) and 4 = 5^/50(5, = 504). In the further numerical 
calculations, we set 6s = 10"^ and 6c = 5^/50, if there is no other statement. 



In FiglHl the dependence of the Born cross sections, full one-loop 0{aeiu) corrected 
cross sections and the corresponding relative corrections of processes pp/pp — > 
X1X2 + at the LHC and the Tevatron on the gaugino mass parameter M2 are 
depicted. There we take the input parameters as m^o=300 GeV, Msusy=350 GeV, 
/i=550 GeV and ^4^=450 GeV, with tan/3 = 4, tan/5 = 15 and tan/5 = 40, respec- 
tively. We can see that the Born and electroweak corrected cross sections decrease 
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Figure 10: The dependence of the full one-loop electroweak corrected cross sections 
forpp/pp — > XiX2+^ processes at the Tevatron and the LHC as the functions of the 
cutoff 6s with 6c = 6s/50(see Fig[T^a),(c)) and 6c with 6s = 506c(see Fig[T^b),(d)), 
respectively. 
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Figure 11: The dependence of the Born, full one-loop electroweak corrected cross 
sections (shown in Fia MlV a.c)] and the corresponding relative corrections 5 (shown 
in FiaYnVh.d)) for the processes pp/pp xf X2 + ^ LHC and the Tevatron on 

the gaugino mass parameter M2. There we take the input parameters as mAo=300 
GeV, Msusy=500 GeV, fi=400 GeV and Af=450 GeV, with tan/3 = A, tan/3 = 15 
and tan /3 = 40 respectively. 
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rapidly to small values with M2 going from 120 GeV to 260 GeV at the Tevatron 
(shown in Fig JTlT c)) and the LHC (shown in Fig llir a)) due to kinematical effects, be- 
cause the masses of the final-state chargino and neutralino are roughly proportional 
to M2. While the relative correction 5(shown in Fig^^b)(d)) increases clearly. At 
the Tevatron the relative correction 6 is in the range from -8.8% to 1.2% with our 
chosen parameters, while at the LHC it can reach -7.9%. Here we have the masses 
of chargino and neutralino in the ranges of m-+ G [111.756 GeV^, 252.504 GeT^], 
m^o e [112.067 GeV, 252.538 GeV]. 

In Figs lT^ a.c). we show the dependence of the Born and the full one-loop elec- 
troweak corrected cross sections of processes pp/pp — > xfx2 + on the SUSY 
soft breaking mass parameter Mg^sy at the LHC and the Tevatron, on the condi- 
tions of m^o=300 GeV, M2=127 GeV, /i=278 GeV and Aj=450 GeV, with tan/3 = 
4(m-+ = 104.428 GeV, m^o = 106.632 GeV), tan/? = 20(m-+=113.353 GeV, m^o = 113.780 
GeV) and tan/5 = 40(m-+=114.577 GeV, m^o=114.839 GeV), respectively. The 
solid curves represent the Born cross sections and the dotted curves represent the 
one-loop electroweak corrected cross sections. We can see the total cross sections 
have the values from 0.34 ph to 1.26 ph and from 3.33 ph to 11.65 ph at the Tevatron 
and the LHC respectively, as the SUSY soft breaking mass parameter M^usy runs 
from 250 GeV to 950 GeV. Figs lT^ b.d) present the one-loop electroweak relative 
corrections as the functions of Mgusy, and we can see the relative corrections gener- 
ally decrease with the increment of Mg^sy And the relative corrections are between 
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Figure 12: The dependence of the Born and the electroweak corrected cross sections 
of process pp/pp — > xfxt + ^ (see in FiaYTWa)fc)) and the corresponding relative 
correctionssee in Fia \l^ b}fd} on the SUSY soft breaking mass parameter M^usy 
the LHC and the Tevatron on the conditions of mAa=300GeV, M2=127 GeV, fi=278 
GeV and Af= 450GeV, with tan/3 = 4, tan/3 = 20 and tan/? = 40, respectively. 
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-9.33% and -7.64% at the Tevatron and between -7.46% and -5.03% at the LHC 
with our chosen parameters. 

We also calculate the Born cross sections, the one-loop electroweak corrected 
cross sections and the corresponding relative corrections as the functions of // (with 
niAo = 300 GeV, Ma = 127 GeV, M^usy = 250 GeV, and Af = 450 GeV, tan/3 = 
4 or 15 and // e [250 GeV, 1000 GeV]), and tan/3 (with tuao = 300 GeV, M2 = 
200 GeV, Msnsy = 350 GeV, Af = 450 GeV, = 550 GeV and tan/3 e [5,40]). 
We find that the results do not depend much on those parameters. The relative 
corrections have the typical values of about —8% and —6% at the Tevatron and the 
LHC, respectively. 

V Summary 

In this paper, we present the calculations of the NLO QCD and the full one-loop 
electroweak corrections to the processes pp/pp — > X1X2 + at the Tevatron and 
the LHC in the framework of the MSSM. In the calculations of both the NLO QCD 
and one-loop electroweak corrections we apply the algorithm of the phase-space slic- 
ing(PSS) method. We analyze the numerical results and investigate the dependence 
of the cross sections and corresponding relative corrections for the processes on sev- 
eral supersymmetric parameters. We find that the NLO QCD corrections generally 
increase the Born cross section, while the electroweak corrections decrease the Born 
cross sections in most of the chosen parameter space. The contributions from the 
NLO QCD corrections make the theoretical predictions nearly independent of the 
renormalization and factorization scales. Our results show that the NLO QCD and 
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electroweak relative corrections typically have the values of about 32%(42%) and 
-8%(-6%) at the Tevatron(LHC), respectively We conclude that the NLO QCD 
and complete one- loop electroweak corrections to the processes pp/pp ^ X1X2 + 
are generally significant and should be considered in high precision analysis. 

Acknowledgments: 

This work was supported in part by the National Natural Science Foundation of 
China and a special fund sponsored by Chinese Academy of Science. 

References 

[1] S. Nandi and X. Tata, Phys. Lett. B129 (1983) 451; A.H. Chamseddine, P. 
Nath and R. Arnowitt, Phys. Lett. B129, (1983) 445; H. Baer and X. Tata, 
Phys. Lett. B155, (1985) 278; H. Baer, K. Hagiwara and X. Tata, Phys. Rev. 
D35, (1987) 1598; P. Nath and R. Arnowitt, Mod. Phys. Lett. A2, (1987) 
331. 

[2] H. Baer and X. Tata, Phys. Rev. D47 (1993) 2739; H. Baer, C. Kao and X. 
Tata, Phys, Rev. D48 (1993) 5175; H. Baer, C.-H. Chen, C. Kao and X. Tata, 
Phys. Rev. D52 (1995) 1565; H. Baer, C.-H. Chen, F. Paige and X. Tata, 
Phys. Rev. D54 (1996) 5866; H. Baer, C.-H. Chen, M. Drees, F. Paige and X. 
Tata, Phys. Rev. D58 (1998) 075008. 

[3] DO Collaboration, S. Abachi et al., Phys. Rev. Lett. 76 (1996) 2228; CDF 
Collaboration, F. Abe et al., Phys. Rev. Lett. 80, (1998) 5275. 



42 



[4] V. Barger and Chung Kao, Phys. Rev. D60 (1999) 115015, Fermilab-Pub- 
98/342-T, MADPH-98-1085, hep-ph/9811489l 

[5] W. Beenakker, R. H5pker, Nucl.Phys.Proc.Suppl. 51C (1996) 261-266, W. 
Beenakker, R. H"opker, M. Spira, P.M. Zerwas, Nucl. Phys. B492 (1997) 51- 
103, W. Beenakker, M. Kramer, T. Plehn, M. Spira, P.M. Zerwas, NucLPhys. 
B515 (1998) 3-14. 

[6] T. Plehn, Ph.D, Thesis, University Hamburg 1998, DESY-THESIS-1998-24, 
|arXiv:h^-ph/9809319l 

[7] E. L. Berger, T. M. P. Tait, M. Klasen, ANL-HEP-CP-01-002, 
|arXiv:hep-ph/0101164 

[8] W. Beenakker, M. Klasen, M. Kramer, T. Plehn, M. Spira, P.M. Zerwas, Phys. 
Rev. Lett. 83 (1999) 3780-3783. 

[9] S. Catani and M.H. Seymour, Nucl. Phys. B485(1997) 291 and Erratum ibid 
B510(197)503. 

[10] K. Fabricius, I. Schmitt, G. Kramer and G. Schierholz, Z. Phys. Cll, 
315(1981); G. Kramer and B. Lampe, Fortschr. Phys. 37, 161(1989); W. T. 
Giele and E.W.N. Glover, Phys. Rev. D46, 1980(1992); W. T. Giele, E.W.N. 
Glover and D.A. Kosower, Nucl. Phys. B403, 633 (1993). 

[11] A. Denner, H. Eck, O. Hahn and J. Kiiblbeck, Nucl. Phys. B387 (1992)467. 

[12] G. P. Lepage, J. Comput. Phys. 27 (1978) 192. 

43 



[13] G. Altarelli and G. Parisi, Nucl. Phus. B126 298 (1997). 

[14] B.W. Harris, J.F. Owens, Phys. Rev. D65(2002)094032. 

[15] Thomas Hahn and Christian Schappacher, Comput. Phys. Commun. 
143(2002) 54, arXiv:hep-ph/0105349l 

[16] J. F. Gunion, H. E. Haber, NucL Phys. B272 (1986) 1. 

[17] J. Pumphn, D.R. Stump, J. Huston, H.L. Lai, P. Nadolsky, W.K. Tung, JHEP 
0207 (2002) 012. 

[18] J. Kiiblbeck, M. Bohm, A. Denner, Comput. Phys. Commun. 60 (1990) 165; 
T. Hahn, Comput. Phys. Commun. 140 (2001) 418; T. Hahn, C. Schappacher, 
Comput. Phys. Commun. 143 (2002) 54; T. Hahn, M. Perez- Victoria, Comput. 
Phys. Commun. 118 (1999) 153. 

[19] G. J. Oldenborgh, Comput. Phys. Commun. 66 (1991) 1; T. Hahn, Acta Phys. 
Polon. B30 (1999) 3469. 

[20] D. M. Copper, D. R. T. Jones, and P. van Nieuwenhuizen, NucL Phys. B167, 
479 (1980); W. Siegel, Phys. Lett. B84, (1979) 193. 

[21] D. A. Ross and J. C. Taylor, NucL Phys. B51, (1979)25. 

[22] A. Denner, Fortschr. Phys. 41 (1993) 307 . 

[23] H. Eberl, M. Kincel, W. Majerotto, Y. Yamada, Phys. Rev. D64(2001) 115013, 
|arXiv:hep-ph/0104109l 

44 



[24] A. Denner, T. Sack, Nucl. Phys. B347, (1990) 203. 

[25] B. A. Kniehl, A. Pilaftsis, Nucl. Phys. B474 (1996) 286. 

[26] J. Guasch, W. Hollik, J. Sola, JHEP 10 (2002) 040, |arXiv:hep-ph/0207364| 



W. Hollik, E. Kraus, M. Roth, C. Rupp, K. Sibold and D. Stockinger, Nucl. 



Phys. B639 (2002) 3, |arXiv:hep-ph/0204350 



[27] D. Pierce and A. Papadopoulos, Phys. Rev. D47 (1993)222; Zhang Ren- You, 
Ma Wen-Gan, Wan Lang-Hui and Jiang Yi, Phys. Rev. D65 (2002)075018. 

[28] A. Sirlin, Phys. Rev. D22(1980) 971; W.J. Marciano and A. Sirlin, Phys. Rev. 
D22(1980) 2695; A. Sirlin and W.J. Marciano, Nucl. Phys. B189(1981) 442. 

[29] C. Weber, H. Eberl, W. Majerotto, Phys. Lett. B572 (2003) 56, 
|arXiv:hi^^/0 305250| 

H. Eberl, M. Kincel, W. Majerotto and Y. Yamada, Nucl. Phys. B625 (2002) 



372, |arXiv:hep-ph/011130"3 



[31] K. Kovarfk, C. Weber, H. Eberl, W. Majerotto, Phys.Lett. B591 (2004) 242- 
254. 

[32] A. de Riijula, R. Petronzio, and A. Savoy-Navarro, Nucl. Phys. B154 (1979) 
394. 

[33] Particle Data Group, Eur. Phys. J. C15 2000. 
[34] S. Eidelman, et al., Phys. Lett. B592 (2004) 1. 

45 



[35] A.D. Martin, R.G. Roberts, W.J. Stirling, and R.S. Thorne, Eur. Phys. J. C39 
(2005) 155, |arXiv:hep-ph/0411Q40 



[36] S. Heinemeyer, W. HoUik, G. Weiglein, Phys. Lett. B455 (1999) 179. 

[37] H. Haber, G. Kane, Phys. Rep. 117 (1985) 75; J. Gunion, H. Haber, Nud. 
Phys. B272 (1986) 1. 



46 



